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Abstract 

We show that Green's functions of second order differential operators 
with singular or unbounded coefficients can have an anomalous behaviour 
in comparison to the well-known properties of the Green's functions of 
operators with bounded coefficients. We discuss some consequences of 
, such an anomalous short or long distance behaviour for a diffusion and 

■ wave propagation in an inhomogeneous medium. 

1 Introduction 

' , A wave propagation with a constant speed is an approximation to the real sit- 

^ nation when the wave propagates in a medium with its characteristics varying 

, in space. A similar approximation is apphed when considering a diffusion. In 

general, the speed of the diffusion can vary in space. If this variation is slow 
then one could believe that its effects are negligible. It has been observed some 
^ ' time ago ( see [1] and references quoted there) that a diffusion in strongly in- 

homogeneous materials can be anomalous. This happens in particular with a 
heat convection in a turbulent medium [2] [3]. In general, it is rather difficult 
to investigate these problems because exact solutions are not available and ap- 
proximations to equations with varying coefficients are not reliable. We discuss 
here the proper time method for a representation of the Green's functions. We 
represent solutions of the equation 

AGe = 26 (1) 

for the Green's function of an elliptic operator A in terms of its heat kernel. 
In this way the Green's function is expressed by a diffusion. A solution of 
such an equation allows to determine the wave propagation if we consider equa- 
tion (l)either as an analog of the Hclmholtz equation for the propagation of 
monochromatic waves (then all the coordinates in (1) are spatial) or continue 



1 



analytically the solution in one coordinate (identified with time) to imagi- 
nary values (we denote such a continuation by Gf)- Then, the Green's function 
G_F is known as the Fcynman propagator. Its real part is a sum of the conven- 
tional advanced and retarded propagators. In general, the Feynman propagator 
is not relevant to classical field theory. However, quantized electromagnetic field 
is necessary for a description of photons moving in an active medium which can 
be described as a dielectric or magnetic material [4] [5] [6] . 

If the coefficients of A are bounded regular functions then it is known [7] that 
there is a minor effect of varying coefficients on the local behaviour of the 
Green's functions. We show that if the coefficients can grow or are singular then 
they have a profound effect on the Green's functions. Then, the change of the 
behaviour of the Green's functions has impact on propagation of disturbances 
in the medium. 

Let us mention some applications of our results to the wave equations in 

physics. Mechanical waves can be derived from Euler and continuity equations 
[8] . For example, the wave equation for the pressure p reads 

ud^p = Vp-^Vp (2) 

where i> is the compressibility and p is the density . An inhomogeneous equa- 
tion whose solution is expressed by the Green's function describes the effect of 
external forces acting upon the medium [9]. There is an equation similar to 
eq.(2) for the scalar potential (j) defined by the velocity v = —Vcj). 

Next, we consider electromagnetic waves. There is a remarkable similarity 
between wave equations in electrodynamics and in fluid dynamics. It follows 
from the Maxwell equations that [10] 

e^t^E = Vfi-^ X V X E - a* J (3) 

where is magnetic susceptibility and e is the dielectric permittivity. The 
equation for the magnetic field can be written in the form 

fid^U = -Ve"^ X V X H V X e'^J (4) 

We could also introduce scalar (p and vector A potentials as 

B = V X A 

E = -V(/) - dtA 

Then, in the static case it follows from Maxwell equations that the scalar po- 
tential (I) satisfies the equation 

djedjcp = —p (5) 

where p is the charge density. The vector potential A is a solution of the wave 
equation similar to eq.(4) [5] [6]. 
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The operator on the rhs of eqs.(3)ancl (4) is mixing the components of the 
electric and magnetic fields. It takes a simpler form if e, ^ and J do not depend 
on one coordinate xd and we are looking for solutions {Ed or Hd) which do 
not depend on xd- Then, cq.(3) reads 

ed^Eo - dj^-'djEn = -dtJo (6) 

Taking the Fourier transform in time of eq.(6) we obtain the Helmholtz equation 

-ew^E^ - djU'^djEf) = -iujJ^ (7) 

The quantum version of this model describes a quantum electromagnetic field 
interacting with a classical source. The correlation functions of the photon field 
and the scattering matrix for photons scattered on the source (or produced by 
the source) are described by the Feynman propagator. Eqs.(6)-(7) are reduced 
to a two-dimensional space. In the next section we discuss the scalar wave equa- 
tion (2) in D space-time dimensions {D — 1 = d space dimension) for greater 
generality. The scalar wave equation is often considered as a good approxima- 
tion to the vector one. The Poisson-typc equation (5) can be considered as 
a special case of eq.(6) when the time derivative is neglected. In this paper 
we discuss quite unusual (unphysical) models for the space-dependent coeffi- 
cients in eqs.(2)-(7). However, our results show that the conventional wisdom 
about spreading of forces in hydrodynamics and electrodynamics may require a 
modification in strongly inhomogeneous media. 



2 Green's functions of second order elliptic dif- 
ferential operators 

In application to the models (2)- (7) we discuss the following equation for the 
Green's functions 

{Bd^ + djAdj)G = 2S (8) 

where can be either a space coordinate, an imaginary time or the real time. In 
the latter case B should be negative. Let us first discuss a special case B = 1. 
Then, 

(d^ + djAdj)GE = 2AGe = 26 (9) 
As an auxiliary tool for the wave equation we consider the diffusion equation 

d.p = \djAdjP (10) 

We consider the fundamental solution Pr of eq.(lO) (the transition function 
of a Markov process). The operator A defined on the Ihs of eq.(9) is a non- 
negative symmetric operator in the Hilbert space of square integrable functions. 
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Hence, it has the unique self-adjoint extension. We can define exp(— r^) and 
subsequently as an integral kernel of 
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dr exp(— T^) 



(the equality holds true on the domain where makes sense; concerning 

Hilbcrt space methods for Green's fimctions see [11]). ,4 is a sum of two com- 
muting operators. Hence, the kernel of exp{—TA) is a product of kernels of 
these operators. In this way the solution of eq.(9) can be expressed by Pr 

Ge{xo,:x.;xq,x.') = j^^duj /g^ (irP(r, x,x') exp(-ia;^T)exp(iw(a;o - x'^)) 
= /o°°<^'r^('r>x,x')(27rT)-5 exp(-^(a;o - x'^f) 

(11) 

Let us denote by qs(x) the diffusion process defined by the transition function 
P (eq.(lO)). Then, the Green's function (8) can be expressed in the form 

GB(a;o,x; a;[,,x') = J du J dr 

E[exp - Ld^ B{c[s{x))dsj exp{iu}{xo - a;o))5(qT(x) - x')] ^^'^^ 

where E[..] denotes the expectation value over the diffusion process qx ([13][14]) 
defined by the transition function P of eq.(lO) . We cannot find a general 
solution of eqs.(8)-(10) but restrict ourselves to 

A{x) = fc|x|" (13) 

where is a constant. 

We can compute the transition function P(t, 0,y) for a diffusion starting 
from a point x = 0. Then, the solution for P in D space-time dimensions is 

P{T,0,y)=K{kT)-'^eM-^\yf) (14) 
where K is a normalization constant resulting from J dyP = 1 and 

13 = 2 -a (15) 

J = {D - 1){2 - a)-' (16) 

The requirement of a normalization of P imposes the condition a < 2. 
We can calculate the variance of the diffusion path 

E[qx(0)2] = j dxP(r,0,x)x2 = Cit^ (17) 
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where Ci is a constant (which can be calculated and depends on the constants 
entering eq.(14), it is finite if a < 2). Hence, for a > we have a superdiffusive 
behaviour whereas for a < a subdiffusivc one. There is no contradiction 
with the diffusion equation (7) which determines the diffusion coefficient as 
limr^oE[T~^q^] because A(x) ^ as x — > in the superdiffusive case and 
^(x) — ^ oo as X ^ in the subdiffusivc behaviour. In eq.(ll) we can calculate 
directly the r integral if 7 + ^ > 1. We obtain 

Gsixo, 0; x'o, x) = C^^Eixo, 0; x'^, x)-t+5 (18) 
with a certain constant C2 (depending on the parameters of eq.(14)) where 

7-i = (D-2 + |)(2-a)-i (19) 

and 

nE{xo,0;x'o,^) = {xo-x'of + ^(1 - ^T^l^f (20) 

The function in the integral (11) is intcgrablc at small r for any "/ ii Qe > 
and for large t if 7 + ^ > 1, i.e. D — 2 > In the hmiting case of 7 = ^ the 
integral (11) is logarithmically divergent for large r. The divergence means that 
the operator A in eq.(9) is not positive definite. We can still define the Green's 
function (9) in this limiting case (as a distribution) subtracting the divergent 
constant from the integral (11) (or imposing the condition J fdx = on the 
test functions /). Then, the Green's function is proportional to Infi^;- 
The special case of eq.(9) is worth mentioning 

djAdjOs = 26 (21) 

The solution is 

Ge(x.,x') =/o°°drP(T,x,x') (22) 

The integral (22) is convergent if 7 > 1 {D — 1 > 2 — a). Inserting the explicit 
solution (14) we obtain 

G£;(0,x) =C|xl-^+3-« (23) 

In the case D — 3 = —a again the logarithm will appear on the rhs of eq.(23). 

The limiting case a ^ 2 is also interesting because in this limit the formulas 
(14)-(20) lose their meaning. In order to study what happens for a = 2 we have 
found the diffusion process whose sample paths give the transition function P 
(solving eq.(lO)). The square of it has the form (qr(x))^ ~ x^ exp(ciT + C2&t) 
where br is the Brownian motion. It follows that qT(0) = is the only solution 
starting from (there is the zero solution also for < a < 2 but in addition a 
non-trivial one (14)). This result can explain why the formulae (14)-(18) have 
no limit as a — > 2. Now, instead of the power-law behaviour (17) we obtain 
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-E[(qr(x))^] = x^exp(cT)( where c > is a certain constant). We suggest that 
when a — > 2 then the power hke relation between q and r resulting from eqs.(14) 
and (17) is replaced by an exponential one following from an exponential increase 
of (such a relation follows also from an explicit calculation of the transition 
function for the radial part of the process qr). 

3 The wave propagation and other apphcations 

We can continue eq.(20) into the real time . The result determines the Feynman 
propagator [12] 

Gpit, 0; t', x) = C2{^F + ie)-^+^ (24) 

where 

np = -{t-t'f + \{i-'^)-M^ (25) 

Let us note some characteristic features of this behaviour. For x = we 
obtain 

Grit, 0; t', 0) = Cs\t - t'|-2(o-2+t)(2-a)-i (26) 
Gpit, 0; t' , 0) can have arbitrarily high singularity (when a ^ 2). lit = t' then 
GF(t,x;t,0) = C4|x|-°+2-t (27) 

Hence, the Green's function is integrable in x. We can see that for a > 

the Green's function is more singular than the one with constant coefficients, 
whereas for a < it is less singular. 

The Feynman propagator does not appear in classical theory of wave propa- 
gation (except of the Feynman- Wheeler electrodynamics) because it propagates 
the waves from the past and from the future as its real part is a sum of the 
advanced and retarded propagators 

^Gf = \{Ga + Gr) (28) 

However, any solution G of eq.(8) {xq = it) is a sum G = Gp + U oiGp and a 
solution U of the wave equation 

{Bd^, - d^Adi)U = 

Choosing a regular U with some boundary conditions (in particular at spatial 
infinity) we can obtain a Green's function with the same singularity and appro- 
priate boundary conditions. 

The real part of Gp is a fractional derivative od 5{Q,f)- If we send a signal 

from at time t then its arrival at time t' at the point x can be calculated 
from the equation Qp{t,0;t' ,yi) =0. Let us note that the wave front does not 
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propagate with a constant speed but the velocity depends on time (or space). 
So, we obtain from Clp 



=2p-'aHl-^)^\t-t'\^ (29) 



dt 

For the more general equation (8) 

(30) 



GE{xo,yi;x'o,-x.') = / dr£'[exp - (xq - 4)2(4 /J^ 5(q^(x))ds) ^ 
(47r/;B(q,(x))rf5)-i^(q.(x)-x')] 



The case a = in eq.(13) {A = k) is also interesting. Then, qs(x) = x + b^, 
where b denotes the Brownian motion [13] [14] (in such a case we know P(r, x, x') 
explicitly). 

The expectation value E[.] in eq.(30) cannot be calculated exactly. We have 
to resort to some approximations. If S is a bounded slowly varying frmction 
then approximately 

S(q,(x))2.i?(x)^i?(0) (31) 
In such a case we obtain the formula (20) for Ge{xo, 0; Xg, x') with 

J^b(xo,0;4,x') = B{Q)-\xo-x'of + ^(1 - (32) 

If B like A has a powerlike behaviour 

B(x) = ^Ixl"" (33) 

then the approximation (31) cannot be justified. Let us note that from eq.(14) 
it follows that 

P(T,0,y) = P(l,0,T-^y) (34) 

Hence, 

q3(0)^T^^(0) (35) 

Inserting this approximate equality (this equality is exact for the Brownian 
motion b)into eq.(30) and making the approximation 

/ B{ci,{0))ds ^ E[ [ B{ciM]ds = C 
Jo Jo 

we obtain 

Ge{xo,0;x'„,^) = KC-i j;;^ drr--^-'^-^ 
exp(-^r-^(.o-.^)^-AM!) (36) 

We cannot compute the integral (36) exactly. Let us consider some special cases. 
If xo = Xq then 

GB(0,0;0,x)~|x|-^+2-t-f (37) 
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If X = x' 

Ge{xo,0;x'q,0)c^\xo-x'o\-'' (38) 

where 

i/ = (2(£)-2) + a + cr)(2-a + a)-^ (39) 

The case a = and a > could be considered as a reahstic approximation 
to eq.(2) describing the density of a fluid concentrated at x = ( then p(x) 
is decreasing from the origin as |x|^"). In electrodynamics a scale invariant 
dielectric permittivity is unrealistic. We could consider 

e = eo + 6|x|« (40) 

where eo is a constant. In application to eq.(5) a source p produces the electric 
potential 

0(y) = -liGEPKy) (41) 

Hence, for a static point source located at x = it follows from eq.(23) that 

Hy) = Cslyf-"-" (42) 

is valid for small distances if a < and for large distances if a > (Z) = 4 
in our world; eq.(42) modifies the conventional ^ Coulomb law). We obtain a 
similar formula from eq.(4) for a magnetic field of a monochromatic wave. 

The function B in eq.(8) has a meaning of the magnetic susceptibility for 
eq.(4). Hence, if 

fi^^io + b\xf (43) 

where is a constant then B{qs{0))ds ~ /xqt + ct^^'^ . In such a case the 
approximations made at eq.(30) hold true for small distances (resulting from 
small r) if a is negative and large distances if a is positive. The same argument 
holds true when applied to eq.(3) with the dielectric permittivity of the form 
(40). Then, the results of sec. 2 concerning scale invariant coefficients apply to 
e of the form (40) if we consider small distances in the case a < and large 
distances if a > 0. 

Finally, we consider QED in an inhomogeneous medium. The quantum 
electromagnetic field is defined by its time-ordered vacuum correlation functions 
which arc determined by the Feynman propagator Gp. As an example of its 
applications wc consider a scattering of photons on an external source J. The 
conventional formula for the S'-matrix can be applied (see [15]) 

.S = exp(- / JGFJ)exp{— i J A(-)j)exp(-j j A^+h) (44) 

where A'-~^ is the part of the vector potential linear in the creation opera- 
tors, whereas is the annihilation part (positive energy solution of the wave 
equation). In particular, if the source is of the form J{t, x) = f{t)S{x) then the 
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probability that the source wiU cause no photon emission ( a preservation of the 
initial vacuum [15]) is 

P(m) = exp(- J duj\f{u)\^GF{u)) (45) 

where tilde denotes the Fourier transform and Gpi^o) is the Fourier transform 
of the Green's function (38). If f{t) = cos{u)Qt) then 

In P(m) ~ uj^~'^d{0) ~ Tw^"^ (46) 

where 6{0) is interpreted as the duration Tof the signal J (this interpretation 
follows from the integral / JGpJ in cq.(44) if T is large in comparison to the 
period 2^ ). We can sec that the inhomogcncity of the medium could be mea- 
sured by a photon counting experiments showing whether ui''~^ is linear in ui or 
not. 



4 Summary 

The potential of a point charge is determined by the Poisson equation (5) with 

p as a (5 source. The electric field can be obtained from eq.(3). The spatial 
singularity and the decrease at infinity of the potential and of the electric field 
are determined by the behaviour of the Green's functions. We have shown that 
this behaviour is not universal and may depend on the properties of the medium 
in which the charge is embedded. Analogous formulae are known for acoustics 
and fluid dynamics expressing the intensity of the waves produced by a given 
source. The laws governing the production of such waves will be modified if 
the fiuid density or its compressibility is a growing or singular function of x. In 
our examples there is always a distinguished point in Ge(x, x') identified as an 
origin x = of the coordinate system. This point can be distinguished by the 
location of the source. The special role of the origin will disappear if coefficients 
in the differential operators are random variables invariant under translations ( 
we discussed such random Green's functions in [16]). 
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